Diffusive persistence and the 'sign-time' distribution 
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We present a new method for extracting the persistence exponent 9 for the diffusion equation, 
based on the distribution P of 'sign-times'. With the aid of a numerically verified Ansatz for P we 
derive an exact formula for 9 in arbitrary spatial dimension d. Our results are in excellent agreement 
with previous numerical studies. Furthermore, our results indicate a qualitative change in P above 
d ~ 36, signalling the existence of a sharp change in the ergodic properties of the diffusion field. 



PACS numbers: 05.40.+J, 82.20.-w 

In the past few years there has been much interest in 
calculating the persistence properties for a wide range of 
simple model systems. Examples are the diffusion equa- 
tion 0]], the Ising model with Glauber dynamics 0j and 
its Potts model generalization ||, interface kinetics jjj, 
phase ordering pi, and the voter model ||. Perhaps the 
simplest and most generic system is the first - the dif- 
fusion equation. Naively one might expect that every- 
thing is known about such a classical system. However, 
one needs only reflect upon its intimate relation to the 
rich Burgers model of turbulence (obtained via a sim- 
ple non-linear transformation) to appreciate the potential 
complexity of diffusion physics. 

This complexity was again uncovered by studies of dif- 
fusive persistence [Q. The persistence exponent 6 for 
this case is defined as follows. Consider the determinis- 
tic diffusion equation evolving a random initial condition 
(usually created from an uncorrelated gaussian distribu- 
tion). Then consider the probability q(t) that the diffu- 
sion field at a given site has never changed sign. One 
finds numerically that this probability decays with time 
in a power-law fashion, with an exponent 8, whose value 
is not a simple rational number. There is no analytic pre- 
diction for 8 with the exception of the results from the 
'independent interval approximation' (IIA) which are in 
good agreement with numerical work in spatial dimen- 
sion d — 1 , but fare less well in higher dimensions |Q . It 
is easy to find applications for diffusive persistence due 
to the ubiquitous presence of diffusion physics. Exam- 
ples include: survival of reactants in reaction kinetics, 
and more general survival probabilities in systems with 
a field slaved to a diffusion process. 

In this Letter we shall present an exact analytic form 
for 8. The key to our derivation is that one may ob- 
tain 8 by studying a more general quantity; namely the 
distribution P of 'sign-times' r(t) (to be defined below), 
which has also been recently introduced by Dornic and 
Godreche [§1 . This distribution may be shown to have an 
exact scaling form for all t: P dr = f(r/t)(dr/t). Nu- 
merically the scaling function / is found to be extremely 



simple. Using this form of / as an Ansatz, allows an exact 
determination of 8. Our prediction deviates by 3% from 
the numerically determined value in d = 1, but lies well 
within the error bars of the simulation results in d = 2 
and d = 3. We shall discuss the reasons for the slight 
d = 1 discrepancy toward the end of the Letter. We also 
find that the qualitative nature of P goes through a sharp 
transition at d c ~ 36. 

We consider the evolution of a scalar field <j)(r, t) which 
satisfies 



d t 4> = d\j< 



(i) 



with initial condition <j>(r, 0) = i/j(r), where the field tjj is 
an uncorrelated random variable described by a gaussian 
distribution 



-(1/2A) / d d r V>(r) 



(2) 



(3) 



R[i>] ~ exp 

The solution of Eq. (0) has the form 

0(r,t) = J d d r' g(r-r',t) ^(r') 

where g(r, t) = {AnDty 1 / 2 exp(-r 2 /4L><) is the heat 
kernel. 

The most obvious way to approach the persistence 
problem is to directly calculate the probability that the 
field at a given site (r = say) has never changed sign. 
This amounts to the evaluation of 



?(*)=(n^(o,o) 



(4) 



where 8{z) is the Heaviside step function |J. Apart 
from the IIA (which is difficult to systematically improve) 
there does not seem to be any possibility of calculating 
this average. 

Let us now focus our attention on 'sign-times' r(t) de- 
fined as 
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t 

r(t) = J dt'9(0(O,t')) , (5) 
o 

so that r(t)/t is the fraction of time t in which the field 
at r = was positive. Note that this is a much easier 
object to handle than q(t) since we do not require the 
positivity of <f> for a continuous succession of times. We 
simply follow the evolution of <j> and record how often it is 
positive. Now, the sign-time r(t) is actually a functional 
of the random variable ip, and as such is described by 
some probability distribution P[t, t] (see also Ref. ||). 
We can immediately list a few properties of P. First, 
it is defined in the domain [0, t]. Second, it is symmet- 
ric about t = t/2 (since it does not matter whether we 
study how often <j) is positive, or how often <p is nega- 
tive). Third, and most important, the behaviour of P 
near r = or t directly furnishes us with the exponent 
9. This is clear, since the probability for r to be in the 
vicinity of either or t is nothing more than q(t) defined 
above. 

We refer the reader to Appendix A in which we prove 
an important fourth property of P. Namely, that it as- 
sumes the scaling form P[r, t] dr — /(r/i)(dr/i) for any 
non-zero time t. [Note that P does not depend on the 
model parameters D and A.] Given the symmetry of P 
about t = t/2, we may rewrite the scaling form as 

P[r,t] dr = g(x(l - x))dx , (6) 

where x = r/t and we have assumed / to be analytic 
around x = 1/2. The third property listed above imposes 
that g(y) ~ y B ~ x for y <C 1, where y = x(l — x). So we 
have reduced the persistence problem to that of calculat- 
ing the tail of the 'sign-time' distribution. This is still a 
formidable task, as it involves calculating arbitrarily high 
moments of the distribution (see Appendix A). As an al- 
ternative strategy we make an Ansatz. Namely, that the 
small-y form for g actually holds for all y in the available 
range y S [0,1/2]. This very simple Ansatz was both 
suggested and confirmed to us by the results of our nu- 
merical work, which we now briefly describe. (Note also 
that this Ansatz appears naturally within the framework 
of the IIA §.) 

Following previous work 0, we model the diffusion 
equation by a discrete space-time process 

&(i + 1) = &(*)+ a -&(*)] , (7) 
j 

where the sum is over nearest neighbours of i on a d- 
dimensional hypercubic lattice. The parameter a is cho- 
sen to be 1 /(2c?) . The initial value of each <fii(0) is drawn 
from a gaussian distribution. Simulations are performed 
on large lattices (N ~ 2 20 sites) for times up to t ~ 2 10 
with several independent runs. We measure two quanti- 
ties during the simulation. First, we record the number 



n(t) of those sites at which the field has never changed 
sign. The ratio of n/N ~ t (due to the self-averaging 
of the system) . We also record the sign-time for each site 
and thus construct the histograms for P[t, t]. Our re- 
sults for the former quantity are in agreement with those 
of Ref. p|. In Figs. 1 and 2 we show the histograms for 
d = 1 and d = 2 respectively. The results are plotted on 
a log-log scale, since the data and the Ansatz are indis- 
tinguishable on a linear scale. One sees that the scaling 
function g indeed varies as a simple power law in y over 
the range y € [0, 1/2]. The only region in which there is 
a deviation from this power-law behaviour is near y = 0. 
This is purely a (time)-lattice effect, since the histogram 
has a finite number of bins given by the number of time- 
steps. As t increases, more and more bins are supplied 
near x = (where the function is integrably singular) 
and the expected power-law behaviour stretches closer 
and closer to the origin. We stress that the power-law 
behaviour near x — is guaranteed by the existence of 
the persistence exponent 8. Our a priori knowledge of g 
is limited in precisely that region of x (G [e, 1/2]) where 
the numerical results show a clear single power-law be- 
haviour. We have added straight-lines on the plots, which 
have slopes of 6 — 1, with values of 9 taken from Eq.(|Tl|) 
below. 
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FIG. 1. A log- log plot of the sign-time distribution P 
against y for d = 1. The data are taken at times t = 2™ 
with n — 6, 7, 8, 9 from top to bottom. The straight line is a 
power law with slope 6 — 1, 8 = 0.1253. 
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FIG. 2. As above, with d = 2 and 9 = 0.1879. 

We are now in a position to calculate 9. The numerical 
results clearly support P[t, t] = (c/t)[x(l — x)] 9 ^ 1 . This 
distribution contains only two parameters: an amplitude 
c (to be set by normalization), and the exponent 9. The 
latter can in principle be set by the calculation of any 
even moment of P (since the odd moments contain no 
new information, as the odd cumulants are zero). The 
simplest to consider is obviously the second moment \i 2 . 
(Unfortunately it is extremely difficult to calculate any 
even moments above the second.) The integrals over P 
are simply given in terms of the Beta function JIo| . From 
the normalization one fixes 1/e = B(9,9). The calcula- 
tion of the second moment yields 



B( 



2,i 



B 



(1 + 0) 
2(1 + 261) 



(8) 



We may now independently calculate [i 2 from the orig- 
inal definition of the sign-times. Explicitly we have: 

l l 

M 2 = ({r(t)/t) 2 ) R = J dai J da 2 C 2 ( ai t,a 2 t) , (9) 
o o 

where C 2 = (0(<j>(O, aii))0(</>(O, a 2 t))) R . This latter 
quantity may be calculated exactly to give (see Appendix 
for a hint) 
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(10) 



The integrals in (J9J) are easily performed, leaving one with 
the expression 



M2 = j(2-/3) 



P 



2(1-/3) 



where 



2n J (1 + a) 
o 



1 + a 
2aV2 



- 1 



-1/2 



(11) 



(12) 



where ip(z) is the digamma function pcfl , showing that 
the persistence exponent is not a simple rational number, 
but is transcendental. We refer the reader to Table 1 
where values of 9(d) are listed, along with the numerical 
and IIA estimates of Ref. |l]] . 

It is of interest to calculate the large-d form for 9. Then 
one finds: 



where 



/ J dir 





[log(l 



1/2 



11/2 



•O(l) 



= 3.0005 • 



(15) 



(16) 



giving a value of 9(d) ~ (0.1850 •• -)Vd. This is to 
be compared to the result from the IIA which gives 
9iiA(d) ~ (0.1454 • • -)Vd. As can be seen from this large- 
d limit, and also from Table I, the IIA consistently un- 
derestimates the value of 9. 

The fact that 9 passes through unity is very interesting 
as it has a direct consequence for the sign-time distribu- 
tion P. For 9 < 1, P has integrably divergent tails at 
r(t) — > and r(t) — > t. Also, the mean of the P (which 
is at r = t/2) is the least likely value of r. One can 
understand this by considering a given point being lo- 
cated in the centre of a very large positive domain. A 
long time must pass before a negative domain sweeps 
through, thus halting the sign-time clock. However, if 
9 > 1, the distribution becomes convex, and the tails go 
to zero at the end-points of the range [0,i\. Thus the 
mean value r = t/2 is in this case the most likely value 
of t. The dynamic mixing of positive and negative do- 
mains is in this case very efficient. (We shall return to 
this point in our conclusions.) It is therefore important 
to know at what dimension 9 passes through unity. This 
is equivalent to insisting that /3(d c ) — 2/3. Numerical 
evaluation of the integral in Eq. (|l^) using Mathematica, 
yields the value d c = 35.967- • -. The fact that this enor- 
mous dimension plays a physical role in diffusion physics 
is extraordinary at first sight: such is the complexity of 
diffusive persistence. 



This integral may be performed explicitly in one and two 
dimensions, with the results: 

m = ^ (11/8) + ^ (9/8) " W8) " v ' (5/8) ] 



8V2/37T - 1 = 0.20042 • ■ 



(13) 
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Oua 


^sim 
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0.1253- • • 


0.1203- • • 


0.1207(5) 


2 


0.1879- • • 


0.1862- • • 


0.1875(10) 


3 


0.2390 ■ • ■ 


0.2358- •• 


0.2380(15) 


> 1 


0.1850 •••Vd 


0.1454- --Vd 





and 



/3(2) = i[>(5/4)-^(3/4)] 

= 4/71- - 1 = 0.27323- 



(14) 



Table 1: The calculated value of 9 from Eqs.(^Tj) and 
(|l2|), along with the IIA and simulation estimates from 
Ref. ||. 

Before concluding, we wish to make some comments 
regarding the role of the space-time lattice. It is quite 
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possible to set up a calculation of the persistence expo- 
nent, or the sign-time distribution, with the discrete al- 
gorithm (fjj) as a starting point. This formulation has the 
advantage that the direct calculation of q(i), as denned 
by Eq.(0), is well-defined. [In the continuum a micro- 
scopic time cut-off must be introduced to make sense of 
the time slices.] We have also pursued the discrete for- 
mulation in an attempt to calculate P[r,t]. As with the 
continuum case, it is extremely difficult to calculate any 
even moment above the second. However, we find the in- 
teresting result that the second moment contains strong 
corrections to scaling; i.e. (t(<) 2 ) ~ t 2 + 0(t). In the 
continuum, each moment has a clean scaling (r(t) n ) ~ t n 
for any non-zero time. The implication of this result, is 
that one should not expect true scaling from numerical 
work (based on (J?])) until very late times. This is espe- 
cially true in low dimensions where the lattice Laplacian 
is much weaker than its continuum counterpart. This ef- 
fect is already apparent in the histograms shown in Figs. 
1 and 2. One sees that the tail (for small y) bends away 
from its power-law form due to the discrete sampling on 
a finite time grid. As more time steps are used, this de- 
viation from scaling is pushed to smaller values of y. It 
should be stressed that the direct measurement of per- 
sistence from the ratio n{t)/N is less reliable since the 
measurement is equivalent to sampling P at y = 
precisely in the region most affected by finite time-step 
effects. It is for these reasons that we believe the d = 1 
measurement of 9, as given in Ref. Q, to be the least 
solid. In principle one can avoid these finite time-grid 
effects by sampling data from the exact solution of the 
diffusion equation, as given in (^|) , where t is a real, con- 
tinuous quantity. Such a numerical study would require 
orders of magnitude more computer time than the previ- 
ous studies, but could potentially give a definitive answer 
for the case d = 1. 

In conclusion we have provided an exact form for 
the diffusive persistence exponent 9, with the aid of an 
Ansatz based on numerical observations of the sign-time 
distribution P. Naturally, this result is in no way rigor- 
ous. We consider a proof of our Ansatz to be a reachable 
goal, although such a proof involves very technical ma- 
nipulations in n-dimensional geometry [ pT|Jl2[ . We have 
also demonstrated the wider significance of the persis- 
tence exponent in parameterizing the sign-time distribu- 
tion. We believe that this distribution will become an 
important new tool in other persistence- type problems. 
The tails of the distribution contain standard persistence 
information, whilst the body of P gives important infor- 
mation regarding the mixing efficiency of the problem at 
hand. If P is concave (convex), than the mean value is 
the least (most) likely. Such information is important 
when ergodic properties of a system are under investi- 
gation. The fact that the curvature of P for the simple 
diffusion equation 'flips' at d — d c ~ 36 indicates that 
the dynamics in the phase space of few-body systems 



may be extremely sensitive to exactly how many degrees 
of freedom are considered. 
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APPENDIX A: 

In this appendix we sketch a brief proof of the asser- 
tion that P satisfies the exact scaling form P[T,t]dr = 
f(r/t)(dT/t) for all t. First, we note that the sign-time 
distribution can be written as P[r,t] = (S(t — t^))r, 
where is the implicit function of ip given in (||). Thus, 



P[r,t] 



J (exp(-to«v)) H . (Al) 



The average on the rhs of the above expression may be 
re-expressed as a power series in terms of u> n (t^) r ; . The 
n th moment is an n-fold integral over the average of n 
step-functions. Each step-function may be represented 
by an integral, yielding (with an implicit limit of ej, — > 0) 



n 

.1 — 1 n 



The average is easily performed over the gaussian distri- 
bution R[ip). On rcscaling the integration variables we 
find 



R> — - 1 n m 



where M l<m = [2a\ /2 aj, 2 / {ai + a m )] d ' 2 . We see that 
the n th moment scales exactly as t n , which enables us to 
combine w and t as a simple product. Returning to flAl]), 
we may scale u> by t to obtain the desired result. 
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